function thehjk(term,type)
load hermitdata;

if(strcmp(type,'polychaos'))
for jj=1:term
    for kk=1:term
        hjk{jj}{kk}= @(x) hermitijk(2,jj,kk)+x-x;
    end
end

elseif(strcmp(type,'determ'))
for jj=1:term
    for kk=1:term
        hjk{jj}{kk}= @(x) 0+x-x;
    end
end
elseif(strcmp(type,'naive'))
% hjk{1}{1}=@(x) 1+x-x;   hjk{2}{1}=@(x) 1+x-x;  hjk{3}{1}=@(x) 1+x-x;   hjk{4}{1}=@(x) 1+x-x;
%    
% hjk{1}{2}=@(x) 1+x-x;   hjk{2}{2}=@(x) 2+x-x;  hjk{3}{2}=@(x) 2-x+x;   hjk{4}{2}=@(x) 1+x-x;
% 
% hjk{1}{3}=@(x) 2+x-x;   hjk{2}{3}=@(x) 1+x-x;  hjk{3}{3}=@(x) 1+x-x;   hjk{4}{3}=@(x) 2+x-x;
% 
% hjk{1}{4}=@(x) 2+x-x;   hjk{2}{4}=@(x) 1+x-x;  hjk{3}{4}=@(x) 2+x-x;   hjk{4}{4}=@(x) 3+x-x;

hjk{1}{1}=@(x) 3+x-x;   hjk{2}{1}=@(x) 2+x-x;  hjk{3}{1}=@(x) 1+x-x;   hjk{4}{1}=@(x) 4+x-x;
   
hjk{1}{2}=@(x) 2+x-x;   hjk{2}{2}=@(x) 1+x-x;  hjk{3}{2}=@(x) 4-x+x;   hjk{4}{2}=@(x) 3+x-x;

hjk{1}{3}=@(x) 1+x-x;   hjk{2}{3}=@(x) 4+x-x;  hjk{3}{3}=@(x) 3+x-x;   hjk{4}{3}=@(x) 2+x-x;

hjk{1}{4}=@(x) 4+x-x;   hjk{2}{4}=@(x) 3+x-x;  hjk{3}{4}=@(x) 2+x-x;   hjk{4}{4}=@(x) 1+x-x;

end


fname='funchjk';
save(fname,'hjk')